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Abstract 

The infinite-genus limit of the KdV- Whitham equations is derived. The limit in- 



m 
o 

volves special scaling for the associated spectral surface such that the integrated den- 
CN ■ sity of states remains finite as N — > oo (the thermodynamic type limit). The limiting 

integro-differential system describes slow evolution of the density of states and can be 
| regarded as the kinetic equation for soliton gas. 

Keywords: modulation equations, finite-gap potentials, rotation number, thermody- 
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1 Introduction 

As is well known 0, [2j, the modulation (Whitham) equations for integrable systems describe 
gradual variations of spectral Riemann surfaces, whose local structure define x, t- oscillating 
1 finite-gap potentials. The wave dynamics problems leading to the modulation equations are 
typically connected with the study of the long-time evolution of large-scale initial data jlj 
and involve the semiclassical asymptotics of the inverse scattering transform. In this asymp- 
totics, the typical scale of modulations 1/e is necessarily proportional to a global number of 
degrees of freedom in the problem g ^> 1 (we assume no small parameters in the original 
equation, instead, we consider large-scale ~ 1/e initial data that can be approximated by g- 
soliton or g-gap potential). Then, asymptotically, the solution manifests itself as a modu- 
lated A-gap potential, where A <C g • The whole point of the traditional modulation 
theory, thus, is to reduce a complicated system with many degrees of freedom to a simpler 
one with a few degrees of freedom. 



The modulation equations, however, can be considered on their own outside any relation 
to the semiclassical limit in the initial value problem. Then, at least formally, the question 
of possibility of the infinite-genus limit for the modulation equations can be risen. In this 
paper, I show that this limit, along with nontrivial mathematical meaning, has natural 
physical interpretation. 

The paper is dedicated to the study of the limit as A — > oo for the A-gap averaged 
Korteweg - de Vries (KdV) equation. The KdV equation is taken in the form 

u t - 6uu x + u xxx = . (1) 

The key in the limiting transition for the modulation equations is the scaling of the spectral 
surface introduced by Venakides in the study of the continuum limit of theta-functions 
[7j. The scaling implies that the essential part of the spectrum of the Shrodinger operator 
associated with lies in a finite interval , say, (—1,0) and the genus of the problem A ^> 1. 



Now, with the width of gaps being 0(1/N) and the bands being exponentially narrow, 
0(exp(— N)), the limit as N — > oo for the KdV-Whitham equations is shown to exist (Section 
5). The existence of the limit for the Whitham equations is ultimately connected with the 
existence of the limit for the quantity called a 'rotation number' (or 'integrated density of 
states') introduced for almost periodic potentials by Johnson and Moser in [S]. Since the 
established limit preserves fmiteness of the density of states as N — > oo, it can be called the 
thermodynamic limit. The limit of this type has been considered recently in [0] in connection 
with stochastic description of infinite-gap potentials. 

The main result of the current paper is presented in Section 6 where the thermodynamic 
limit for the KdV- Whitham equations is derived in the form of a nonlinear transport equation 
()5ip for the density of states. The connection between 'density' and 'velocity' in (joTj) is 
given by the integral equation (|52jl . I suggest that the obtained limit describes macroscopic 
dynamics of the 'thermodynamic limit of stochastic soliton lattices' considered in [0j and can 
be regarded as the kinetic equation for solitons. This interpretation is also supported by the 
exact agreement of the small density asymptotics of the resulting system with Zakharov's 
kinetic equation for a rarefied soliton gas [TUj . 



2 The Whitham equations: basic definitions 

In this section, we give a brief account on some basic definitions and notations from the 
finite-gap theory that will be used in this paper. The KdV- Whitham system of genus N can 
be represented as one generating equation in the Flashka-Forest-McLaughlin form [2j 

d T dp N = d x dq N , (2) 

X = ex , T = et , e< 1, 

where dp^ and dq^ are the meromorphic differentials on the hyperelliptic Riemann surface 
of genus iV 

2N+1 

T : R\E) = J] ( E ~ E i) > E E C (3) 

Ei < E 2 < . . . < E 2 n < E 2 n+i , Ej G R 

with cuts along bands [E 2 j~i, E 2 j\. The canonical basis of cycles is organized as follows: the 
ctj-cycle surrounds the j'-th cut clockwise on the upper sheet, and the (3j- cycle is canonically 
conjugated to a/s such that the closed contour (3j starts at E 2 j , goes to +oo on the upper 
sheet and returns to E 2 j on the lower sheet . 

Then the meromorphic differentials are defiened by: 

d PN (E) = + N ~ 1 R(E) ° dE, fdp N (E) = 0, j = l,...,iV (4) 

Pi 

rpN+l I „ rpN j_ , „ i 2N+1 
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fdq N (E) = 0, j = l,...,N. 

The differentials dp^ and dq^ are often referred to as quasimomentum and quasienergy dif- 
ferentials. The integrals of dp^ and dq^ over the a - cycles are known to give the components 
of the wave number and the frequency vectors j2], JI] 

<j> dp N {E) = kj , I dq N {E) = Wj , j = 1,...,N (6) 

The finite-gap potentials associated with the Riemann surface (J3J are known to be the 
multiphase solutions of the KdV equation ([U) 12 , [T3|, [2] 

u N (x,t;X,T)=u N (y 1 ,...,y N ), (7) 

where the real phases 

y. = kjX + Ujt + fj , j = l,...,N. (8) 

are the linear function of 'fast' x and t while the dependence on 'slow' variables X and T, 
enters the potential via the branch points Ej (X, T) of the Riemann surface. The constant 
fj(mod2ir) is the fixed initial phase. The function ujv(y) is 27r-periodic in each phase yj 

Un{vi, ■■■,Vj + 2vr, . . . , y N ) = u N (y 1 , y N ) , (9) 

that is u^(x,t) for fixed Ei,...E 2 n+i is iV-quasiperiodic both in x and t (hereafter we 
consider only incommensurable kjS and oUjS). 

The Whitham equations (J2J) thus describe slow evolution (modulations) of rapidly oscil- 
lating finite-gap potentials. 

The algebro-geometric structure of the finite-gap potentials is given by the basis holo- 
morphic differentials 

N-l j^k 
fc=0 ^(^) 

and the functions a,jk(Ex, • • • , E2N+1) are determined by the normalization over the a-cycles 

*Pj = Sjk, (11) 



"A- 



while the integrals over the /3-cycles give the entries of the Riemann matrix 

S y =/Vi- (12) 



The wavenumbers fcj and the frequences cjj are expressed in terms of by the formulas 

k = -47ri J B" 1 a JV _i , (13) 

2N+1 

u = (ui, . . . ,u) N ) = -87riS _1 (ajv-i E ^' + 2a iV-2) • 
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3 Rotation numbers and the density of states 

We consider the Schrodinger equation 

(-% x + q(x))i/> = E,l>, x£R. (14) 

with an almost periodic potential q(x). The potential q(x) has an important characteristics, 
the rotation number, which was defined by Johnson and Moser [S] as 

a(E) = lim -axgM/(x,E) + ii/;(x,E)) , E e R (15) 

x — >oo % 

As shown in jS] , this agrees with the value ttM where M{E) is the integrated density of states 
defined by 

M{E)= lim ^iMl, (16) 
|6-a|-»oo b - a 

where u(a, b; E) is the number of eigenvalues Ej < E in the Dirichlet problem on a < x < b: 
ip(a,E) = ip(b,E) = 0. Thus, 

M{E) = —a(E) . (17) 

71 

If, as a particular case, q(x) is iV-gap potential , q(x) = un(x), such that — 1 < un < then 
the rotation number coincides with the real part of quasimomentum [Hj: 



a N (E) =Re J dp N (E') , E e (-1, 0) 



Then, with the aid of © we have 



I M{E) 

a N (E) = (19) 

X M(E) E 

- ^ ij+ / dp N (E') if Eeba,nd M , 
J— E2M-1 

where M(E) is the number of the band nearest to E from the left; M(E) < N , M(0) = 
N , M(-l) = 0. 

We note that the formula (fT5|l agrees with the general statement in [H] that the values 
of the function 2a(E) if E G {gap} belong to the frequency- module of the almost periodic 
potential. One can also see that the function a'(E) is positive and is supported on the 
spectrum of the finite-gap potential. 

Since the finite-gap potential ((7j), (jHJ) is an almost periodic function in t as well as in x 
one can also formally introduce a temporal analog of the rotation number 

E 

(3 N (E) =ReJ dq N (E') , (20) 
-1 
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where dq^{E) is the quasienergy differential (J3J). 

We also introduce the full integrated density of states k which is a number: 

1 N 

n=Af N (0) = —J2kj. (21) 

Z7F 3=1 

One can see that the full integrated density of states has for the finite-gap potentials the 
natural significance of the mean number of waves per unit length. 
Analogously, we introduce 

1 1 N 

O = -fo(0) = — 5>,- . (22) 
2tt . =1 



4 Continuum limit of the spectrum and its thermody- 
namic nature 

Let the nontrivial (finite band) part of the spectrum lies in the interval (—1, 0) of the real E- 
axis: ajv(— 1) = 0, «Ar(0) = ttk. For N ^> 1 we consider the special band/gap distribution 
(scaling): 

gaps(E) ~ bands(E) ~ exp (—AT) N > 1 , Be(-1,0). (23) 

Following Venakides we introduce the lattice of points 

1 sa 77! > 772 > • • • > r/jv ~ , (24) 

where 

2 



are centers of bands. 

We define two positive continuous functions on [0, 1]: 

1. The normalized density of bands <p{jj): 

, . , number of lattice points in (77, 77 + drj) 
(p{r])dr) « — . 

That is, 

^ fe)= N(r li -r li+1 ) +0 ^ ) - /^)*»= 1 - I* = "^(0,1). 

2. The normalized logarithmic band width 7(77): 

7(77;) = ~ log 5, + O(l) , 5, = £ 2j - ^ . (26) 

The functions ^p(r]) and 7(17) asymptotically define the local structure of the Riemann surface 
T (jSJ) for ^> 1. In other words, instead of 2 A" + 1 discrete parameters Ej we have two 
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continuous functions of rj on (—1,0), which do not depend on x,t on Ax ~ At ~ 1. On a 
larger scale, Ax ~ At ~ 1/e , e < 1, in the spirit of the modulation theory, we put 



(p = <p(r);X,T) 7 = 7(r?;X,T) 



(27) 



where X = ex, T = et. These functions are now subject to the evolution equations which 
naturally should be derived as the continuum limit of the Whitham equations (J2J) considered 
for the same exponential spectral scaling (|23|) . 

In the next section we will show that the scaling (}2~3*j) implies the following asymptotic 
behaviour for the wavenumbers and the frequences (these asymptotics can also be derived 
directly from (j§J)): 



kj ~ Uj ~ 



log 5, 



<4> 



iV> 1, 



:28i 



I dp N (E')~ J dq N (E')<0(^), E E (E 2 j_i, E 2 j) , (29) 

which provides finiteness of the rotation number ajv(-E') (and of the integrated density of 
states A/at(S)) as N — > oo: 



a(E) = lim 

N— >oc 



M(E) 



lim 



< oo . 



M < N . 



(30) 



Due to this property we call the scaling (|23|) and the corresponding limit as N — > oo the 
thermodynamic ones. This type of limits in the finite-gap theory has been introduced in Pj. 
Hereafter we will denote the limit of some functional J[un] as N — > oo on the thermodynamic 
spectral scaling as 

T- lim f[u N ] 



It should be noted that the thermodynamic spectral scaling (|23j) appears when one considers 
the semiclassical asymptotics of the spectrum for periodic potentials ^3], @], However, in 
the works just cited, N ~ 1/e, which is the standard relation in the asymptotics of this type 
(see also j3]). Contrastingly, in the thermodynamic limiting transition, the large parameter 
N defining the scaling is contained only in the band-gap structure and is not connected with 
the parameter 1/e, which scales modulations of the finite-gap potential UN(x,t). 



5 The thermodynamic limit for the rotation numbers: 
basic integral equations 

Our aim now is to deduce the thermodynamic limit for the modulation equation (j2J). For 
that, we need to evaluate the thermodynamic limits of the meromorphic differentials dp^ 
and dqjsf. First we make an estimate for their imaginary parts on the thermodynamic scaling 
(El: 
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E G any band 



Im J dp N (E') 



-1 



(31) 



dp N (E') 



E 2j 



<0( 



' N 



E e gap. 



Therefore, the imaginary part of dp^(E) vanishes as iV — > oo and we have for E e (—1, 0): 



1 M(E) 

T-\imp N (E) = a (E) = T-lim- £ k 3 - 

2 3=1 



Analogously, we get 



M (E) 



T-\imq N {E) = (3(E) = T-lim- ]T 



UJ 



3 i 



(32) 



(33) 



i=i 



where p(E) = T-lim p N (E). 

Now we evaluate the wave numbers kj and the frequencies cjj for the chosen spectral 
scaling. For that, we will make use of the following Venakides observation [7j: for the 
spectral distribution (|23|) the following approximation for the basis holomorphic differentials 
is valid: 



*l>3 



N 

n 



2tt 



27V+1 



1/2 



dE. 



(34) 



n 



One can see that 



\ k=l 

(E - E 2 j-i) (E - Ey) = (E + 7$ + 0(e 



Then an extensive cancellation becomes possible in (J34"|) so that the hyperelliptic integrals 
in (jH$, ()12|) turn into elementary ones. As a result, we have the following asymptotics (see 
|7| for details): 



*7» - 



^ + ?7j 



1-3 ' 



N. 



0-j,N-l 



2tt 



(35) 
(36) 



2iV+l 

a 3,N-l >, -Ej + 2aj,N-2 
3=1 

Now we represent the general relationships (fT3*|) between the real and imaginary periods of 
theta function in the form convenient for the thermodynamic limiting transition 



£>k = — 47riajv_! 



Buj 



2N+1 

^7Ti(a iV _i E 3 + 2a iV-2) • 



(37) 
(38) 



It follows from the asymptotics (|35|). (|3l))l that to provide a balance of terms in (J37j) . (|55Jl the 
following scaling is necessary : fc, ~ Uj ~ 0(1/N). Therefore, we introduce two continuous 
on (0, 1) functions k{rj) and uj{rj) by 



(39) 



N J N 

Then, substituting (jHSJ), (JHSJ) , (JSHJ), into (pT7j) . (jHHjl we get two algebraic systems for k(r)i 



AT j 
JV x 

— losr 



Vi - Vj 



Vi + Vj 
Vi - Vj 



KVj) +l(Vi)HVi 



-27lT]i . 



u(Vj)+l(Vi)u(r)i) = 8?7i 



(40) 
(41) 



Vi + Vj 

Using ([39)1 we represent the relationships (|32p. (|33j) for the thermodynamic limit of the 
rotation numbers in a continuum form: 

2da(-r] 2 ) = (p(rj)k(r])dr) , 2df3(-rf) = ^(77)^(77)^/7, (42) 

Now one readily obtains the thermodynamic limit of (|4(J|) . (J41j) : 



1 

/log 


77 - p 





V + H 



a'(— \j?}d\i + cr(rj)c£'(— rf 



7i 

2 



1 

/log 


77-// 





77 + /X 



^/3'(-fi 2 )dfi + a(r])(3'(- V ' 
77 



-27T77 2 



where 



(7(77) 



7(5) 



(43) 



(44) 



(45) 



We see that the thermodynamic limit of the rotation numbers is determined by the only 
function a{rj) defined on (0, 1) (instead of the original two: <p{rf) and 7(77)). 

The integral equations (J43|) . (J44j) have been derived in [_0 in connection with the es- 
tablishing the thermodynamic limit for the stochastic processes generated by the finite-gap 
potentials (stochastic soliton lattices). We note that these integral equations also appear in 
[16 where they determine the Lax-Levermore type minimizer for the iV-soliton solution with 
randomly distributed soliton phases as N — > 00. It is clear that there should be a direct 
connection between the 'stochastic' version of the Lax-Levermore variational problem and 
the thermodynamic limit of the rotation number. 



6 Thermodynamic limit of the Whitham equations 

The modulation equation (J2J) describe slow X,T- deformations of the spectral surface (J3J). 
Now we are going to apply the thermodynamic limiting transition to (0). As follows from 
dSH) ; the local thermodynamic limits are 

T-\imdp N (E) = da(E) T-\imdq N (E) = d(3(E) . (46) 
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We put that on a larger scale 



da = da(-r] 2 ;X,T) dp = dp(-r] 2 ; X,T) . 
Then the thermodynamic limit of the Whitham system (j2J) takes the form: 

<9tc/ = dx/3' , 



(47) 



(45 



where the functional dependence j3'[a'(— r] 2 ; X, T)] is given by ()43j) . (03). The thermody- 
namic limit of the Whitham equations, therefore, describes evolution of the density of states 
a'(E)/n. 

Since da(E) > one can introduce a distribution function f(rj;X,T) by 



f( m X,T)d V =-da(-r] 2 ;X,T) 

7T 



± 

/ fdr, 



K . 



Introducing 



one can represent 



s(r,;X,T) 



a'(-^;X,T) 
in the form of a nonlinear transport equation: 



d T f = d x (sf), 

where the integral equation connecting s with / is obtained by excluding a(rj) from 



(49) 



(50) 



(51) 



s(r)) 



V 



1 

/log 


v-n 





7] + fl 



f(fi)[s(ji) - s(rj)]dfx 



(52) 



Here, of course, s(rj) = s(rj; X, T), f(r)) = f(r]] X, T). In the case of small integrated density 
of states /j C 1 we get 



4 
V 



1 

/log 


Tj-fl 








(53) 



which is Zakharov's expression for the velocity of the 'trial' soliton moving through rarefied 
soliton gas (TUJ - The integro-differential system (JoTj) . (jUJ, therefore, can be regarded as the 
finite-density generalisation of the Zakharov kinetic equation for solitons. 

Now we discuss briefly what happens to the underlying finite-gap oscillating microstruc- 
ture under the thermodynamic limiting transition. It has been shown in jU] that the ther- 
modynamic limit of the finite-gap potentials should be treated as a random process. More 
accurately, the thermodynamic limit is defined not for the individual potential un{x) but 
for the stationary random process generated by considering the associated torus equipped 
with the Haar measure. Such processes were called in jTTj stochastic soliton lattices. The 
stochastic soliton lattice v{x\4>) is defined by 



U N (x\(f)) =U N (...,kjX + (f)j, ...), j 



N. 



(54) 
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where = (0i, . . . , 0at) is a random value uniformly distributed on Tor^. It was shown in 
jH] that the thermodynamic limit of the stochastic soliton lattice 

T- lim ujsf(x\<p) = ((x — Xj) (55) 

represents a random process defined on the infinite discrete set of points {xj]j G Z} dis- 
tributed by Poisson on R. Since the band/gap ratio for the chosen spectral scaling is propor- 
tional to Ne~ N , it vanishes as iV — > oo and the thermodynamic limit of the stochastic soliton 
lattices (|53j) can be associated with one- dimensional soliton gas. This also supports our in- 
terpretation of the thermodynamic limit of the Whitham equations as the kinetic equation 
for solitons. 



7 Averaged conservation laws 

Integrating (J48|) by E over the interval (—1,0) we obtain the thermodynamic limit of the 
wave conservation law: 

d T K = d x £l ■ (56) 

where ttk is the full density of waves (integrated wavenumber) (}2~T]) and nVt is the corre- 
sponding 'integrated frequency' (see (j221))- 

The quasimomentum Pn{E) is known to be the generating function for the Kruskal 
integrals I ( k N) jUJ, [Ej: 



oo AN) 
1 1, 



^(^ = 2^(1 + 17^). -£»1, (57) 

1 L 1 L 1 

Iq = lim — / u N (x)dx , l[ N ' = lim — / -u%(x)dx , . . . (58) 

L^oo L J L— >oo L J 2 


Then the thermodynamic limit for the Kruskal integrals is calculated as 

P k (X,T) = T-Um/f » = T-lim (r^^I^-E)^ 

22fc+l " 



l) k+1 J(-E) 2Jl f 1 a'(E;X,T)dE 



-i 

i 



T (-l) fc+1 J V 2k+1 f(m X, T)d V , k = 0, 1, 2, . . . (59) 



(2fc+ 1)tt 

22fc+l 

~ 2k + 

u 

Thus, the infinite set of the independent averaged conservation laws takes the form 

d T P k = d x Q k , k = 0,1,2 ... , (60) 

where the 'fluxes' Qk{X,T) are found analogously to (|59|) from the thermodynamic limit for 
the quasienergy differential Po^l : 

22fc+i \ 

Q k (X, T) = ^—^(-1)^ j r, 2k+1 f( V] X, T)s{ m X, T)d V . (61) 

o 
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In particular, it follows from (|58p. ()59|) that the thermodynamic limits for the first two KdV 
moments are 

i lfi i 

u(X, T) = -2 J 77/(77; X, T)dr/ , ^(X, T) = - J V 3 f( V ; X, T)d V . (62) 


We note in conclusion, that it is clear that the procedure of the thermodynamic limiting 
transition can be extended to other integrable Whitham hierarchies. Another implication 
of the developed theory is the thermodynamic limit for the Krichever algebro-geometrical 
procedure of integration of the Whitham equations jTHj, [HJ. This result will be published 
in a separate paper [T9] . 
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